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We have constructed microscopic, structurally-relaxed atomistic models of Si/Si02 superlattices. 
The structural distortion and oxidation-state characteristics of the interface Si atoms are examined 
in detail. The role played by the interface Si suboxides in raising the band gap and producing 
dispersionless energy bands is established. The suboxide atoms are shown to generate an abrupt 
interface layer about 1.60 A thick. Bandstructure and optical- absorption calculations at the Fermi 
Golden rule level are used to demonstrate that increasing confinement leads to (a) direct bandgaps, 
(b) a blue shift in the spectrum, and (c) an enhancement of the absorption intensity in the threshold- 
energy region. Some aspects of this behaviour appear not only in the symmetry direction associated 
with the superlattice axis, but also in the orthogonal plane directions. We conclude that, in contrast 
to Si/Ge, Si/Si02 superlattices show clear optical enhancement and a shift of the optical spectrum 
into the region useful for many opto-electronic applications. 

PACS numbers: 78.66.Jg, 68.65.+g, 71.23.Cq 



I. INTRODUCTION 

The initial interest in light-emitting Si-based nanos- 
tructures has lead to a number of important experiments 
establishing that Si/Si02 superlaiticfiSiJjSLs) show en- 
hanced, blue-shifted luminescence .BB'cIqijIj While the lu- 
minescence pattern is more complex in these systems 
than in others, the blue shift was found to correlate with 
decreasing Si-layer thickness. This simple relation be- 
tween the silicon-layer thickness and the luminescence 
peak is of great interest for applications such as Si-based 
light-emitting diodes (Si-LEDs). All reported SL energy 
peaks are in the lower part of the visible spectrum — 
the highest reported value being 2.3 eV (540 nm), i.e., 
green apd the lowest 1.2 eV (1030 nm), in the near 
infrared. u These SLs are really multiple Si quantum wells 
(MQW) , the silicon oxide layers playing the role of bar- 
riers. The thickness of the Si quantum wells, igj, is 
the critical parameter. A fixed igi would simply set the 
colour of the Si-LED, while MQWs with a range of Lsi 
would span a range of colours in the luminescence. 

Other systems containing confined |Silicon structures 
(besides MQWs) are, e.g., porous SiEna (consisting of 
quasi one-dimensional structureis),, silicon nanoclusters 
in Si02 matrices,EI nanocrystalst£l or dislocation loops, 
and qjUantum-dot structures made from implantation of 
bororJid or other ions £3 The choice of a particular struc- 
ture, e.g., for Si-LED applications, depends on many fac- 
tors: stability over time, optical efficiency at room tem- 
perature, experimental reproductiveness, facility to ac- 
cept n- or p- type dopants (e.g., for n-p junctions), ease 
of incorporation in ultra-large-scalc-intcgration technol- 
ogy, and production costs. Si/Si02 SLs are stable struc- 
tures, as opposed to porous-Si. In addition, the silicon 
layer thicknesses in SLs are directly related to the energy 



peaks in the Si/Si02 luminescence spectra. This is not 
straighforward in porous silicon. Si clusters or nanocrys- 
tals, where pore dimensions as well as hydrogen concept, 
trations play an uncontrolled role on the energy shift Jj'E^I 

Silicon-based structures have many advantages over 
structures made of other semiconductors: low-cost (as 
compared to any of the IH-V's or H-VTs), non-toxicity, 
practically unlimited availability (in contrast to germa- 
nium) and benefits from decades of experience in pu- 
rification, growth, and device fabrication. However, the 
indirect energy gap (~1.1 eV) in bulk crystalline silicon 
(c-Si) makes it unsuitable for optoelectronic applications. 
Silica (Si02) is another key material of the microelec- 
tronics industry; it has a bandgap of '^9 eV. Optical- 
fiber technologies and metal-oxide-semiconductor field- 
effect transistors (MOSFETs) are based on (high-quafity) 
silica. Molecular-beam epitaxy and chemical vapor depo- 
sition provide the needed growth technology fop^. Si . aj ftd. 
Si02. It is possible to combine crystalline siliconBllJlla't^l 
and Si02 to produce structured materials having chem- 
ically pure, sharp, defect-free interfaces. The degree of 
advancement of the fabrication technology is such that 
the enhanced luminescence in SL/Si02 SLs cannot be ex- 
plained only in terms of defectsQ and/or residual hydro- 
gen atoms filling unsaturated Si dangling lupnds — the 
latter being referred to as Pb-type centeralZl — located 
at the SL interfaces, as in the well-understood enhanced 
luminescence of hydrogenated a-Si. 

The objective of the present article is to consider the 
detailed microscopic structure of the Si02/Si/Si02 dou- 
ble interface structure and provide a first-principles un- 
derstanding of the emission of light from Si/Si02 SLs. 
Only a few atomistic i^DEI-caleulations on Si/Si02 SLs 
have been reportedEa'tBE^EilOS Many more would have 
been available were it not of the (naturally-occuring) 
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amorphous structure of Si02: amorphous structures re- 
quire large supercells to get physically relevant results. 
Another important issue is the need to model the Si/Si02 
interface so as to correctly incorporate the known expaii^ 
imcntal details. Recent core-level shift measurementsEil 
provide-|details of the suboxide (partially oxidized) Si 
atoms.E3 Further, the abruptness of the Si/Si02 inter- 
face has been established from transmission electronic 
microscopy (TEM) experiments; values of thejititerface 
width as low as 5 A have been reportedE'EEl Based 
on these observations, realistic Si/Si02 ii 
els have been designed by several workers. lil 
Suboxide Si atoms in most of these models are distributed 
within three atomic layers of the Si/Si02 interface, corre- 
sponding to the lowest experimental interface thickness. 
Such interface models are needed for first-principles mod- 
elling of MOSFETs, and for SL structures, where multiple 
Si/Si02 interfaces are present. 

An early (and tractable) Si/SijQ2 interface model was 
proposed by Herman and Batra.E3 The large lattice mis- 
match between Si02 and Si was accommodated by set- 
ting the /3-cristobalite Si02 unit cell diagonally on the 
diamond-like c-Si unit cell. An oxygen atom was included 
in the interface to saturate the dangling bonds, resulting 
in a crystalline model of the interface. Another simple 
crystalline model, involving a bridge oxygen at the ifu 
terface, was introduced by Tit and Dharma-wardana.LJ 
These models were studied using a variety of methods 
— tight- binding (TB),E3i23_full-potential, linear- muffin- 
tin-orbitals (FP-LMT0)E3cul and linearized-augmented- 
plane- waves (FP-LAPW)E3 The dispersionless charac- 
ter of the bandstructure in the growth axis has been 
confirmed within all three theoretical approaches, thus 
demonstrating the existence of strong confinement. How- 
ever, the nature of the energy gap — from both LMTO 
and LAPW calculations — is still indirect. Further- 
more, these crystalline models show that the light ab- 
sorption is quite dependent on the details of bonding 
and interface structure (bondlengths, angles and chem- 
ical species), emphasizing ihe need for more realistic, 
structurally-relaxed models.E3 

Several structurally-relaxed models have been investi- 
gated by Kageshima and Shiraishi using first-principles 
methods. Their models were constructed starting from 
/3-cristobalite as well as a-quartz Si02 layers superposed 
onto c-Si layers with different possibilities for the Si/SiOe, 
interface (such as hydrogen atoms at dangling bonds) ;t2l 
the atomic sites were then structurally relaxed. The cal- 
culations indicate that the energy gaps are indeed direct, 
and that interfacial Si-OH bonds are possible candidates 
for the light-emitting enhancement in SLs. However, the 
models are not consistent with the observed suboxide 
Si atomic distributions at the Si(001)/SiO2 interfaces. 
For instance, no Si atoms bonded t£L. three oxygens are 
present. Tit and Dharma-wardanaO have constructed 
a partially-relaxed model (PRM) starting from a struc- 
turally relaxed Si(001)/SiO2 interface strttcture due to 
Pasquarello, Hybertsen, and Car (PHC).EZl This inter- 



face model contains all three suboxide Si atomic species. 
Hydrogen atoms were used by PHC to terminate the 
surface. In the Tit and Dharma-wardana model, the H 
atoms were removed and the Si/Si02 interface structure 
was converted into a Si02/Si/Si02 double interface SL 
structure while preserving local tetrahedral bonding to 
obtain the PRM structure. Within the TB approach, 
the energy gap of the PRM was shown to be direct and 
enhancement of th£J3ptical absorption (as compared to c- 
Si) was confirmed.Ej This is further described in the next 
section, since the fully relaxed models (FRMs) discussed 
in the present work are based on this PRM. 

The results presented here go beyond the TB approach 
and were obtained within the projector-augmented wave 
(PAW) theoretical framework. A brief description of the 
theory is given in section II. The SL models were all 
structurally relaxed, and contain no hydrogen atoms at 
the interfaces.|-The interface suboxide Si atoms observed 
in experimenter arise naturally in the model. Calcula- 
tions were carried out for different Si-layer thicknesses, in 
order to assess the effect of confinement on the electronic 
and optical properties. Supplementary models have been 
constructed to clarify the role of the suboxide atoms on 
the SL optical properties. We first review the theoreti- 
cal methods (section II), then focus on the four central 
issues needed to understand the Si/Si02 luminescence 
properties: construction of realistic interface models (sec- 
tion HI), quantum confinement (section IV), role of the 
suboxide atoms (section V), and optical effects associ- 
ated with increased confinement (section VI) . Our study 
thus provides a complete, microscopic picture of the lu- 
minescence properties of Si/Si02 SLs. 



II. COMPUTATIONAL DETAILS 

The electronic-structure calculations were carried oiit 
using the Vienna ab initio simulation paokage (VASP)Ej 
using the "frozen-core" PAW approacb-ifj The overall 



framework is density- functional theoryEfo (DFT) within 
the local-density approximation (LDA).e2I 

The "frozen-core" PAW is a simplp: form of the general 
PAW method introduced by BlochlO Blochl's method is 
an extension of the usual LAPWc3 approach. Hence, 
the PAW method formally bridges the LAPW to the 
ultrasoft-pseudopotentials (US-PP) in order to combine 
the precision of the former and the rapidity (for larger 
systems) of the latter. The PAW method has another ad- 
vantage over the usual implementation of US-PP, essen- 
tial for optical calculations: it avoids correcting for spa- 
tial nonlocality effects in typical pseudopotentials when 
evaluating the momentum operator p. Further details 
can be found in the article of Adolph et alE3 

Matrix elements of the momentum operator p are 
needed for calculating interband optical effects. We de- 
scribe the main steps that lead to the calculation of the 
absorption coefficient, starting from PAW solutions of the 
Kohn-Sham equations. The PAW approach rests on the 
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following linear transformation: 



This relates the (calculated) pseudo wave function \^n) 
to the (now corrected) all-electrons wave function IvPat). 
The index specifies the atomic site, angular momentum 
numbers and reference energy. The two functions \(I)n) 
and \4>n) are respectively the pseudo, and all-electron 
wave function of a reference atom. They are forced to 
overlap outside a given core region. The functions \pi) 
are the projector functions characteristic of the PAW 
method. Thus, the three functions |0jv), |</>jv) and \pi) 
constitute the frozen-core PAW data, being set prior to 
self-consistent field calculations. The projector functions 
\pi) are constructed so as to remain dual to the pseudo 
wave functions, to fulfill generalized orthogonality con- 
straints and to remain (approximately) complete (see 
Ref. ^ for full definitions). 

The application of the above linear transformation to 
any operator A within the PAW approach has been de- 
scribed by Blochl [cf. eq. (11) of Ref. |4|. If A is the 
momentum operator p (in a certain direction defined by 
the polarization vector), we have 

(*Ar|p|*M> = (*w|p|*A/) + 

The imaginary part of the dielectric function, e^, can be 
determined by using the Fermi Golden pile within the 
Coulomb gauge; the expression becomesiQ 
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where n = 27re/m. The function /„ is the Fermi distribu- 
tion and n\v\^ m) are the PAW matrix elements. The 
whole expression corresponds to the probability per unit 
volume for a transition of an electron from the valence 
band state m) to the conduction band state I^Pm) to 
occur. |— I 

The tetrahedron methodea is used to evaluate ei{E). 
The joint density of states, which determines the inter- 
band transitions 5{Em — Em — E) and the optical matrix 
elements |(4'Ar|p|^M)P, are computed on each tetrahe- 
dron (i.e., 1/4 X the sum of the matrix elements on the 
four corners of each tetrahedron). The real pact is then 
obtained using the Kramers-Kronig relation.Efl Since the 
dielectric function is the square of the complex refractive 
index, (e^ + iu) = (rir + irii)"^ , the absorption coefRcient 
becomes 



with c the speed of light in vacuum, h Planck's constant, 
and E the photon energy. 

Electron-hole (e-h, excitonic) interactions were not in- 
cluded in the calculations, as they would be in, say, solu- 
tions of the Bethc-Salpeter equation. The size of our sys- 
tems prohibits such complete optical calculations, which 
are feasible only for very small systems (a few atoms). 
These additional effects would enhance the results from 
interband transitions since e-h intSLactions generally in- 
crease the absorption at the onset .Ej 



III. CONSTRUCTION OF THE STRUCTURAL 
MODELS 

Recent core-level shift experimentsil have revealed the 
presence of all possible oxidation states for Si atoms 
in Si/Si02 structures such as SLs, that is Si+", where 
n = 0, 1, 2, 3, 4 is the charge found within each Si Wigner- 
Seitz sphere. Si" and Si"*""^ are the charge states of Si 
found in bulk-Si and bulk-Si02. The suboxide (subO) 
Si atoms with n = 1,2,3 are found at the interface. 
Slightly larger distributions for the subO Si"'"'^ densities, 
as compared to those for Si"*"^ and Si"*"^, were reported in 
these experiments. Microscopic Si/Si02 interface models 
should be consistent with experiments in closely repro- 
ducing the density distributions of all subO Si atoms. 

As mentioned above, the Si/Si02 SL model struc- 
tures discussed in the present article are based on the. 
Si(001)/SiO2 interface structures obtained by PHC,lZI 
who used the Car-Parrinello method to relax the models 
to their energy minima. It is important to note that the 
PHC models were not designed for the double interface 
structure found in Si02/Si/Si02 SLs but, rather, for a 
single Si/Si02 interface which terminates into the vac- 
uum; this is done by saturating dangling bonds with H 
atoms. Thus, these models de facto contain the essen- 
tial details of atomic positions and charge states at the 
Si/Si02 interface. j—. 

Tit and Dharma-wardanao have generated a Si/Si02 
SL structure starting from one of the PHC models that 
contains an equal distribution of the three subO atoms, 
in (almost complete) accord with experiment. This SL 
model has been constructed by first operating a mirror 
transformation and then a partial rotation of the Si/Si02 
section of the relaxed interface structure, leading to an 
intermediate Si/Si02[mirror]Si02/Si SL structure. Sec- 
ond, some of the Si layers were inverted in order to 
meet the sp'^-bonding topology. The resulting Si/Si02 
SLs structure, fulhi. described in the article of Tit and 
Dharma-wardana,Ej has the following final configuration. 
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The connection between the symbols in this configura- 
tion and the specific atomic layers in the model is shown 
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in Fig. H The letters A, B, C and D correspond to sil- 
icon atomic layers, while the Os are oxygen layers. The 
primes denote layers that depart from the diamond-like- 
Si crystalline arrangement, and [O] corresponds to the 
layer where the mirror operation has been performed. _B 
and D are the Si layers which have been inverted in or- 
der to satisfy the sp'^-bonding topology. The subO Si+" 
atoms, with n = 1,2,3, are distributed within only two 
Si layers, while the Si° atoms are distributed within five 
atomic layers. The embryonic PHC interface model cor- 
responds roughly to one side of the above configuration, 
starting from the arrow up to the right. 

Of course, this construction induces artificial symme- 
tries in the middle of the Si02 layer (more precisely, upon 
and around the [O] layer); this model is thus in essence 
partially relaxed — the PRM referred to earlier. Signif- 
icant information on the electronic and optical proper- 
ties of this model have been extracted j— within the TB 
approach, by Tit and Dharma-wardana,c3 who obtained 
direct energy gaps as well as dispersionless bandstruc- 
tures. Furthermore, the imaginary part of the dielectric 
function was calculated, and then the absorption coeffi- 
cient was deduced. From this calculation, enhancement 
of absorption as welLas blueshift with confinement have 
been demonstrated.L3 

The next obvious step is to relax the PRM, i.e., deter- 
mine the set of positions which leads, via the Hellman- 
Feynman forces, to the lowest total energy. We have 
used the PAW approach described in the previous sec- 
tion to obtain a first fully- relaxed model (FRMl), which 
contains approximately one unit cell of confined-Si. The 
supercell contains 52 Si and 44 O atoms, and has dimen- 
sions 7.675 X 7.675 x 24.621 A^. The relaxation proce- 
dure has been performed with five k points in the reduced 
Brillouin zone (BZ). The energy cutoff was 25.96 Ryd in 
all calculations. The total energy was found to decrease 
by 30.84 eV (0.32 eV per atom) during relaxation. Figure 
^ shows the bondlength distributions before (PRM) and 
after (FRMl) relaxation; the bondlengths are centered 
around the expected values, viz., ~ 1.61 A for Si"0 and 
2.35 A for Si-Si bonds. The shaded boxes in Fig. ^ are 
the distributions of the interface subO Si atoms, while 
the empty boxes are the total distribution bondlengths, 
including subO Si atoms; the shaded boxes remain rel- 
atively unchanged upon relaxation since both interfaces 
were already at their energy minimum, after PHC. The 
main atomic drift during relaxation occurs in the cen- 
ter of the Si and Si02 layers, i.e., near the [O] and the 
A layer of the configuration discussed above. Interfacial 
Si-0 bondlengths of all subO Si atoms depart from the 
values of Si"*"* in the Si02 layer. The broadening of the 
Si-Si bondlengths is in general much larger than that of 
Si-0; the distortion of the bondlengths are thus mainly 
within the Si layer and at the Si/Si02 interface, i.e., not 
inside the silica layer. This is further discussed below. 
The resulting FRMl is shown in Fig. |l|. 

Additional models having thicker Si wells were con- 
structed in order to examine the role of subO Si layers 



and the effect of confinement on the electronic and opti- 
cal properties. As noted earlier, the FRMl structure con- 
tains approximately one unit cell of confined-Si, i.e., the 
set of layers with charge state Si*^ (bulk-like Si) . By in- 
serting one, then two, additional ABCD Si atomic planes 
(i.e., one Si unit cell, thickness 5.43 A), and relaxing all 
atoms, we generated two additional models - FRM2 and 
FRM3. The FRM2 contains 68 Si atoms while the FRM3 
has 84 Si atoms; both have 44 oxygen atoms, as in the 
FRMl. The total energy variation for the FRM2 dur- 
ing relaxation was found to be only 0.15 eV (i.e., 0.0013 
eV/atom) while for the FRM3, this change is a minuscule 
0.051 eV (i.e., 0.00040 eV/atom). These numbers imply 
that both FRM2 and FRM3 have essentially crystalline 
Si layers. The FRM2 has nine Si*^ atomic planes while 
the FRM3 has thirteen. 

Figure ^ shows the distributions of the Si- 
Si bondlengths in the FRM3, starting from the 
Si(001)/SiO2 interfaces (at the bottom of Fig. |) and 
going towards the centre of the silicon layer along the 
growth axis. The standard deviation (a) from the mean 
value (a;=2.34 A in all atomic layers except at the inter- 
face, where a;=2.29A ) is also given. The diamond-shaped 
symbols correspond to the Si-Si bondlengths for Si° 
atoms while the filled circles are subO Si-Si bondlengths 
at the interfaces. The Si-Si bondlengths depart signif- 
icantly from their crystalline counterparts at the inter- 
faces up to about three atomic layers, where cr=0.019; 
the standard deviation is four times higher at the inter- 
faces than in the fifth atomic layer. This deviation of the 
bondlengths at the Si/Si02 interfaces shows that it is im- 
portant to take relaxation aspects into account. Indeed, 
amorphous Si02 layers in realistic SLs induce strain and 
disorder in the Si layer, as confirmed in Fig. |^. This effect 
could generate localized defects giving rise to efficient ra- 
diative electron-hole recombination. However, the strain 
fields only contribute to the small quasi- momenta regime 
and cannot easily supply the momentum deficit involved 
in the indirect transition of c-Si. Moreover, the relatively 
small size of our supercell models in the x-y directions 
prevents firm conclusions being drawn about the influ- 
ence of this strain on the optical properties. Further 
aspects of the role of interfaces are discussed in section 
V. 

The three subO Si configurations at the interfaces of 
the FRMs are shown in Fig. ^ with their corresponding 
bondlengths. Note that the left and right interfaces in the 
Si02/Si/Si02 SLs are not exactly equivalent; they remain 
independent (during structural relaxation, for instance). 
However, a subO Si on the left interface has a locally 
equivalent subO Si on the right interface, by construction. 
As a consequence, each pair of equivalent subO Si atoms 
have approximately the same bondlengths and angles in 
all the FRMs. As seen in Fig. ^ the bondlengths depart 
from their bulk values, which are 2.35 A for c-Si and 
1.61 A for Si02. In addition, the angles of the subO Si 
tetrahedra vary considerably: the Si-Si-Si angles vary 
from 99° to 125°, the O-Si-Si angles vary from 96° to 
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126°, while all 0-Si-O angles remain around 106°. It is 
thus clear that the subO Si tetrahedra at the interfaces 
are distorted as compared to bulk-Si tetrahedra. 

As discussed later on, the role of subO Si atoms 
was further studied using the following variations of the 
FRM2: First, we removed all Si+^ atoms and attached 
the proper number of hydrogen atoms to neutralize the 
excess charges. The H positions were then relaxed while 
keeping all the silicon and oxygen atoms fixed. This 
structure thus contains Si° atoms and Si"*"" subOs, where 
n = 1,2,3 (i.e., n ^ A). A variation of this structure 
was generated by removing all oxygen atoms and filling 
the Si dangling bonds with H atoms, and again relax- 
ing the H atoms. The final Si-H bondlengths vary from 
1.47 A to 1.53 A, after relaxation. This final structure 
is thus subO-free and contains only Si*^ atoms, except at 
the interface with the vacuum, where hydrogen atoms 
fill the dangling bonds. These three confinement models 
are shown in Fig. they will be referred to as "FRM2" , 
"FRM2/0-H/vacuum" and "FRM2/H/vacuum" , respec- 
tively. 

In Fig. ^we show the BZ of the supercell, the standard 
c-Si diamond BZ, and the high symmetry axes used for 
the bandstructure calculations. We also constructed the 
bulk c-Si structure in a supercell of dimensions similar to 
that of the FRM SLs, so that comparisons can be done 
within the same k space zone scheme. This will be used in 
the next three sections for comparisons of bandstructure 
as well as absorption calculations. 



IV. QUANTUM CONFINEMENT 

In this section, we discuss the nature of the confined 
states in the SLs. We calculated the bandstructures of 
the three SL models — FRMl, FRM2 and FRM3 ~ 
as well as the supplementary FRM2/0-H/vacuum struc- 
ture, cf. Fig. P(b). The latter is the "ultimate" in terms 
of confinement, as the two interfaces with vacuum con- 
stitute infinite potential walls. All bandstructures are 
analysed and compared within equivalent supercell BZ. 
The growth axis of the SLs being the ^-axis, confinement 
effects are expected to take place in the X-R and Z-T 
axis of the BZ (see Fig. ^for axis definitions). 

The bandstructures of the three SLs and their total 
density of states (DOS) are shown in Fig. |^. We find that 
the bandstructures in the growth axis {X-R and Z-T) are 
dispersionless, for all models and all energies. In phys- 
ical terms, dispersionless bandstructures imply infinite 
effective masses, reflecting the strong confinement. The 
DOS have a more abrupt variation in the valence bands 
than in the conduction bands. However, DOS alone are 
not enough to fully understand the optical processes in- 
volving the interband transitions, since the weighting of 
the optical matrix elements is needed. This is further 
discussed in Section VI. 

Let us consider in more details one of the SL mod- 
els, namely the FRM2. We select the R-Z-T-M high- 



symmetry axis where the major features, viz. the relevant 
energy gaps, appear. The bandstructures of the folded 
c-Si structure, the FRM2 SLs structure as well as the 
FRM2/0-H/vacuum structure have been calculated and 
compared. Figure |8|(a) and ||(b) contain a synopsis of 
all calculations for k points along R-Z-T-M . Several 
conclusions can be drawn from these figures. 

By comparing the bands for c-Si, Fig. ||(a), with those 
for the FRM2 SL, Fig. |(b) (sohd lines), we see that fold- 
ing effects cannot by themselves explain the new band- 
structure. The bands in the Z-T directions are totally 
modified by the confinement. Although c-Si always has 
an optically indirect bandgap, this bandgap becomes al- 
most direct in its folded configuration, as can be seen 
from Fig. ^(a), while for the SL the bandgap is un- 
equivocally direct, and significantly increased. Moreover, 
comparing the bands in a direction orthogonal to the 
growth axis, for instance around M in both Fig. ||(a) and 
Fig. ||(b), we see that the valence bands are raised; the 
lowest conduction band in the T-M direction is pushed to 
higher energies. In addition, they exhibit less dispersion 
in the SL than in c-Si, in general. Hence, confinement 
modifies the electronic properties in the growth axis as 
well as in directions orthogonal to it. This is further 
analysed from optical absorption calculations, below. 

We compare in Fig. 0(b) the bandstructures of the 
FRM2 and FRM2/0-H/vacuum models. The positions 
of the Si° atoms, as well as the subO Si in the two struc- 
tures, are identical. The solid lines displays the band- 
structures of the FRM2, while the dots display the bands 
of the FRM2/0-H/ vacuum structure. It is clear that the 
two bandstructures are nearly identical. This calculation 
shows that the Si02 layers act as virtually impenetrable 
barriers. The electronic properties of hypothetical SL 
structures having only subO Si - i.e., no Si^* of Si02 - 
and positioned at the subO Si sites, would give nearly 
identical electronic properties as Si/Si02 SLs, for ener- 
gies close to the band gap. Our calculations show that 
the electronic wavefunctions die out at the suboxide Si 
atoms of the interfaces. Thus the barrier is sharply lo- 
cated just behind the subO Si atoms, and therefore just 
two atomic layers could be used as a barrier without al- 
tering the electronic properties, when energies involved 
(in the device) remain close to the energy gap. The in- 
fiucncc of the interface subO Si atoms is further analysed 
next. 



V. ROLE OF INTERFACES 

In order to assess the role of the subO ions on the elec- 
tronic properties, we calculated the bandstructures of the 
FRM2/H/vacuum structure. Fig. ||(c), which contains no 
subO Si; the dangling bonds have been filled by hydrogen 
atoms and hydrogen atoms (only) have been structurally 
relaxed. Figure ^ summarizes the results. The LDA band 
gap is still direct but significantly lowered, from 0.81 eV 
in the FRM2 to 0.67 eV in the FRM2/H/vacuum struc- 
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ture. Interface reconstruction (as discussed in section 
III and Fig. ||) thus has significant impact on the elec- 
tronic properties. The bands in the plane orthogonal 
to the growth axis are quite different from FRM2; e.g., 
the valence band is lowered near the M and R points, 
becoming similar to the c-Si bandstructure. We thus 
conclude that the subO Si atoms have two effects: (i) 
increase the bandgap and (ii) produce dispersionless va- 
lence bands. The charge states in the subO are respon- 
sible for the increase in the bandgap, while the strongly 
increased valence-band offset leads to essentially disper- 
sionless bands. 



VI. BLUESHIFT AND OPTICAL 
ENHANCEMENT 

The matrix elements entering the calculations of the 
optical properties are often approximated as a constant 
in a given range of energies. However, such an approxi- 
mation is inadequate for elucidating the enhanced lumi- 
nescence in Si/Si02 SLs. This section deals with calcu- 
lating the absorption coefficient within the Fermi Golden 
rule and interband-transition theory. 

Since the Si-layer thickness Lsi is relevant to the 
energy-shift in SLs, this quantity needs first to be de- 
fined. This involves some uncertainty associated with 
the interface thicknesses Lgubo- 

The interface thickness Lsubo was estimated to be 
-1.60 A from our calculations of the subO Si region. 
This is the largest distance (projected onto the z axis) 
between any two subO Si atoms. Hence the upper- 
bound to Lsi are, 11.17 A, 16.58 A and 22.01 A for the 
FRMl, FRM2 and FRM3 respectively, while the lower- 
bound thicknesses are simply Lsi — 2Lsubo- We define 
the Si thickness in the SLs to include the interface subO 
Si atoms as well (corresponding to the upper bound). 
This choice is made since subO Si (Si"*"^, Si"*"^ or Si"'"'^ ) 
atoms contribute to the electronic properties, as do bulk- 
Si atoms (Si° atoms); for instance, we showed above [see 
e.g. Fig. gb)] that the bandstructures of the FRM2 SLs 
and the FRM2/0-H/vacuum systems overlap, and indi- 
cated where the effective barrier begins. 

The band gaps of the FRMl, FRM2 and FRM3 (see 
the bandstructures in Fig. 0) are direct except for the 
FRMl where the band gap is nearly direct with only 
0.12 eV between the direct and indirect transitions. The 
values of the gap are 0.99 eV, 0.81 eV and 0.68 eV for the 
FRMl, FRM2 and FRM3, respectively. The direct tran- 
sition at r for the FRMl equals 1.11 eV. For the FRM2 
and the FRM3, the band gaps are direct and located on 
the whole Z-T axis [see Fig. 0(b) and (c)]. Direct transi- 
tions can thus be achieved between the valence band and 
the conduction band, along the Z — T line of the BZ. We 
thus obtain, under the LDA, a blue shift 

0.68 eV — > 0.81 eV — > 0.99 eV 



with increased confinement 

2.2 nm > 1.7 nm — > 1.1 nm. 

However, these values for the energy gaps are much lower 
than the experimental ones. It is well known that DFT 
within the LDA underestimates the energy gaps of semi- 
conductors and insulators. For c-Si, the DFT-LDA gap 
is approximately 0.6 eV less than the experimental value. 
For the /3-cristobalite phase of Si02, in the group /42(i, 
the DFT-LDA energy gap is 5.8 eV while the experimen- 
tal value is about 3 eV higher, i.e., —9 eV. Approximate, 
but realistic, band gaps can be obtained by adding 0.6 
eV overall: 

1.28 eV — > 1.41 eV — > 1.59 eV. 

(2.2 nm) (1.7 nm) (1.1 nm) 

These energw gaps correspond to the lower bound of the 
experimental and lie in the visible spectrum. However, 
these gap^-Mfl i^itill somewhat lower than the experimen- 
tal values .ua'Q'ffil The discrepancy can be explained by 
recrystallization processes, which lead to the formation 
of nanoclusters that would increase the confinement, and 
correspondingly the measured energy gaps.cZI The analy- 
sis of such a behaviour, which is beyond the scope of the 
present work, would require zero-dimensional-confined 
model structures. 

The absorption of the three SL models and c-Si 
have been calculated both in the diamond-like BZ and 
in the SLs BZ. For c-Si in iha diamond-like BZ, we 
used (20 X 20 X 20) k points,^ while in the SLs BZ, 
(7 X 7 X 2) k points are used. Calculations using more 
k points, viz. (8 x 8 x 3), show that (7 x 7 x 2) is quite 
sufficient to recover the form of the absorption curve for 
all the SL models. The purpose of calculating the ab- 
sorption of c-Si in two different BZ is to estimate errors, 
first due to zone folding effects (which cause round-off 
errors p-Leading to non-absolutely-null transitions at the 
onset)c3 and, second, to the tetrahedron method itself 
which needs large amounts of k points. The broadening 
in the absorption curves has been fixed to 0.015 eV for 
all the absorption curves discussed below, as suggested 
by FugglcEJ 

Figure |l^ shows the absorption results. Panels (a) and 

(b) give an overall view of the absorption curves for the 
z axis in (a) and the x-y plane in (b). Panels (c) and 
(d) show the absorption at the onset, for the z axis in 

(c) and the x-y plane in (d). In all cases, we included 
the absorption of c-Si calculated in the diamond-like BZ, 
as well as the one in the SLs BZ. Direct comparison of 
the two c-Si absorption curves give an estimate of impre- 
cisions due to zone foldings and intermediate number of 
k point effects. It shows that the absorption is slightly 
underestimated in the SL calculations; e.g., in Fig. p^(c), 
the onset of absorption of c-Si in the SLs BZ take place 
at 2.0 eV while in the diamond-like BZ the onset hap- 
pens at the correct value of 2.52 eV (which is for c-Si 
the LDA direct transition at P). This numerical effect 
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cannot be avoided and will arise, as well, in any Si/Si02 
supercell. Hence, all comparison of the SLs absorption 
must be made with c-Si calculated in the equivalent SLs 
BZ, i.e., within equivalent k space zone schemes. 

Since the SLs are fabricated with the objective of 
changing the indirect gap to a direct gap, we now dis- 
cuss the absorption threshold region. Comparison at the 
onset of absorption from Fig. |l^(c) or (d) shows that all 
absorption curves have a lower energy threshold than both 
the c-Si absorption curves, and especially below the one 
calculated in equivalent SLs BZ having equal number of k 
points. That is, the SLs show absorption (and emission) 
in the spectral region above the indirect gap of c-Si and 
below the direct gap of c-Si. This shows that the absorp- 
tion in all confined Si/Si02 SL models is enhanced, com- 
pared to c-Si; the transitions are direct in SLs and have 
an active oscillator strength. For the folded c-Si energy 
bands, the lower bands above the Fermi level, and the 
corresponding oscillator strength, remain dark; in other 
words, the optical matrix elements of the SL BZ of c-Si 
are null up to ~2.0 eV. This result clearly demonstrates 
the enhancement of the absorption (and emission) mech- 
anisms in confined Si structures. Furthermore, upon in- 
spection of the absorption curves in the plane orthogonal 
to the growth axis [Fig. ^0|(d)], we note that the energy 
thresholds of the absorption are all below c-Si: thus, the 
x-y absorption of the FRMs takes place approximately 
at their respective direct energy gaps, and then behave 
in a similar manner, as expected from the similarity of 
the SLs bandstructures, in this plane. 

We examine, finally, the higher-energy region which 
corresponds to the usual direct transition (3-6 eV) in 
c-Si. Even here. Fig. ^(a) demonstrates a blueshift with 
increased confinement in the z axis. The overall ab- 
sorption maxima for FRM1-FRM3 are at 5.28 eV, 4.83 
eV, and 4.71 eV, with intensities of 136, 155 and 162 
(xlO*/cm) respectively. For c-Si in the SLs BZ (to en- 
sure comparable precision in the calculations) the second 
peak, i.e., the maximum, take place at 4.70 eV (with ab- 
sorption equal to 231 x 10'*/cm), while the first peak is 
at 4.0 eV (and with absorption equal to 229 xlO'^/cm). 
We emphazise that there is still a slight blueshift in the 
x-y plane orthogonal to the growth axis, but less pro- 
nounced than in the growth axis. 



VII. CONCLUDING REMARKS 

In this work, the structural, electronic and optical 
properties of Si/Si02 superlattices have been studied on 
the basis of structurally-relaxed models. These SL mod- 
els, which contain no hydrogen atoms at the Si/Si02 in- 
terfaces, exhibit enhanced optical absorption/emission, 
as observed in experiment; this can be attributed to the 
presence of silicon dioxide barriers. In experiments per- 
formed under ultra-high vacuum conditions, the oxidiza- 
tion process would predominantly give rise to oxide bonds 
at the interfaces, but still, few hydrogen atoms are ex- 



pected to be present and fill some of the remaining dan- 
gling bonds. 

Our calculations show that the oxide barriers are cen- 
tral to the optical enhancement in SLs. It is well known 
that hydrogen atoms play a similar role in amorphous 
silicon by filling dangling bonds. This suggests that it 
might also be the case in Si/Si02 SLs, where hydrogen 
atoms fill extra dangling bonds, and thus would amplify 
the optical enhancement effect already exerted by the 
oxide barriers. Further studies are needed to ascertain 
this. 

We have shown that suboxide Si atoms at the interfaces 
act as virtually impenetrable barriers. The active barrier 
thickness thus corresponds to the suboxide Si layer — 
only 1.6 A in our models. The confined Si layer thus 
consists of bulk Si and suboxide Si atoms. Suboxide Si 
atoms at the interfaces modify the electronic properties 
in two manners: they (i) increase the energy gap and (ii) 
lead to dispersionless band structures, which increases 
the transition probabilities. 

Other confinement models (zero-dimensional struc- 
tures) and interface effects will be considered in future 
work. For instance, inclusion of other atomic species — 
such as nitrogen that would generate subnitric Si atoms 
at the Si/Si02 interface — are expected to modify the 
electronic properties, and hence enhance the optical ab- 
sorption/emission spectra. 
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FIG. 1: The unit cell of the fully relaxed SL model (FRMl). 
The configuration of the bulk-like and suboxide Si atomic 
planes is also depicted. The white and black circles are re- 
spectively the positions of Si and O atoms. 
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FIG. 2: Evolution of the Si-O and Si-Si bondlengths from 

(a) the PRM construction by Tit and Dharma-wardana, to 

(b) the fully relaxed structure (FRMl) described in the text. 
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FIG. 3; Si-Si bondlengtii distribution in tlie FRM3, from the 
interface (bottom of figure) towards the center (top of figure) 
of the Si layer. (There are thirteen Si° layers in the FRM3 
and thus six interplanar Si-Si bondlengths starting from both 
interfaces towards the center of the Si layer. 






FIG. 4: Structure of the three suboxide interfacial Si atoms. 
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(a) FRM2; SLs (b) FRM2/0-H/vacuum (c) FRM2/H/vacuum 




FIG. 5: (a) The FRM2 SLs (b) the FRM2/0-H/vacuum 
model (c) the FRM2/H/vacuum model. 




FIG. 6: Definition of the SLs BZ superposed to the diamond- 
like BZ. The principal axis used for bandstructure calculation 
are also depicted. 
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X 20 40 60 80 IOC 

DOS [States/eV] 



FIG. 7: Band structures and density of states (DOS) of the 
three SL models. The DOS are calculated using (7x7x1) 
k-point grid, which corresponds to 144 irreducible tetrahedra. 
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FIG. 8: Band structures of (a) c-Si in the (folded) SL BZ 
(b) comparison between the FRM2 SLs and the FRM2/0- 
H/vacuum bandstructures. 




FIG. 9: Comparison of the bandstructures for 

FRM2/H/vacuum and FRM2 SLs. 
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(a) z axis (b) x-y plane 
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FIG. 10: Absorption coefRcient of the SL models as com- 
pared to c-Si. (a) and (c) are the absorption in the growth 
ajcis; (b) and (d) are the absorption in the plane orthogonal 
to the growth axis. 



